Electrically controlled quantum gates for two-spin qubits in two double quantum dots 
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Exchange-coupled singlet-triplet spin qubits in two gate-defined double quantum dots are considered theo- 
retically. Using charge density operators to describe the double-dot orbital states, we calculate the Coulomb 
couplings between the qubits, and identify optimal bias points for single- and two-qubit operations, as well as 
convenient idle positions. The same intuitive formulation is used to derive dephasing rates of these qubits due 
to the fluctuating charge environment, thereby providing the main considerations for a quantum computation 
architecture that is within current experimental capabilities. 
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Electron spins localized in semiconductor quantum dots 
(QDs) have been vigorously pursued in recent years as poten- 
tial qubit candidates, due to their relative isolation from the 
environment f 1]. The long coherence times measured for sin- 
gle electron spins render their manipulation and readout par- 
ticularly challenging. This has prompted a number of pro- 
posals to encode the logical qubit into two-spin singlet (S) 
and unpolarized triplet (Tq) states [2-4]. GaAs gate-defined 
lateral QDs in particular have served as a useful platform to 
study quantum control and decoherence of exchange-coupled 
spin qubits. The main reason is the high tunability of the two- 
electron exchange interaction, J, that is controlled on a sub- 
nanosecond scale by changing the bias between the two dots. 
In the context of the S—Tq qubit this tunability enables initial- 
ization of the qubit state in a doubly occupied singlet state us- 
ing positive interdot bias, where triplet states are energetically 
inaccessible. By applying a perpendicular magnetic field, the 
two polarized triplet states are Zeeman spUtted, further reduc- 
ing leakages from the logical basis states. Recent experiments 
in gated double dots have implemented spin-echo [5-7] and 
pulse optimization techniques that have pushed two-spin qubit 
dephasing time over 200 fis [7]. 

In this work we consider two proximal S — Tq qubits, each 
formed by two electrons locaUzed in a gate-defined double 
dot, as depicted in Fig. [Ha) f4, 8T. Single qubit Z-rotations are 
fixed by pulsing the interdot bias that controls J, whereas X- 
rotations are generated by a magnetic field gradient across the 
two dots. A promising method to produce such gradient is to 
utilize the hyperfine coupling at the S* — r_|_ degeneracy point 
to dynamically induce unequal nuclear polarization across the 
two dots. Implementing these electrically controlled pump 
cycles, Foletti et al. were recently able to sustain a stable 
Overhauser field gradient of 5h ~ 5/zeV IsHlOfl. The single- 
qubit Hamiltonian can thus be written as T-Lq = B • cr, where 
B = |(5/i, 0, J), and <j is the vector of Pauli spin matrices 
for the pseudospin states S, and Tq. 

When two double dots are placed near each other, their re- 
spective electrons become Coulomb coupled. This coupling 
depends on the orbital states and thus on the spin configura- 
tions of the electrons, resulting in an effective additional ex- 
change interaction for each qubit, and spin dependent entan- 
gling coupling. Below we show that for certain bias points 
the first term turns off the exchange interaction that is oth- 



erwise constantly on. The latter term provides a controlled- 
phase (CPHASE) gate, thus completing a universal set of 
gates that is based solely on accurate control over the inter- 
dot bias at each qubit. In particular, the requirement to con- 
trol the tunnel coupling between the dots, needed to tune J to 
zero i4 lllll . and the need to perform single-qubit Z-rotations 
near the dephasing-susceptible avoided crossing point 11211 . 
are eliminated, and qubit decoupling positions are identified. 
Inter-qubit Coulomb couplings and sweet spots. The 
Coulomb interaction between two qubits, each formed by the 
singlet and triplet spin states of two electrons in a biased dou- 
ble dot, is found by considering the charge distributions of 
the electron orbital states, assuming no inter-qubit tunnel cou- 
pling. Up to a constant term, the interaction Hamiltonian can 
be written in the form 11 1311 : 



Hi: 



-P^al-f5,al-^at 



(1) 



where 



Pt/c = T (VtoTo T ^STo ± VtoS - V'ss) 

1 , 

7 = ^ (-"^ToTo + VsTo + VtoS - Vss) • (2) 

Here the subscript t (c) denotes target (control) qubit, and 
Vij are the Coulomb matrix elements involving the hybridized 
singlet {S) and unpolarized triplet (To) states, where the left 
(right) subscript corresponds the target (control) qubit. For 
nominally identical, equally biased QDs, we have fit = Pc- 
When the inter-qubit distance R is sufficiently large, there is 
no tunneling between the qubits, the two charge distributions 
are separated in space, and the Coulomb interaction between 
them can be described systematically using a multipole expan- 
sion approach f 13, 141. For the same-axis geometry shown in 
Fig.flja), couplings are always positive, and they are plotted in 
Fig- Eb) as a function of the target qubit bias, et, for R = 300 
nm and a fixed control qubit bias, €c, corresponding to the 
anticrossing point. Here and for all results reported in this 
work we have considered B ~ 100 mT {Ez = 2.5/ieV), and 
have modeled the double dots using a quartic potential with 
dot confinement loq = 3meV (a^ « 20 nm), and half interdot 
distance a ~ 2.8aB. While all the Coulomb couplings depend 
on both {et, ec), for the parameter range considered, each of 
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FIG. 1: (color online) (a) Geometry of the two double-dot qubits. We 
denote the left qubit as the target, and the right qubit as the control; 
(b) Two-qubit Coulomb couplings, and exchange interaction vs. tar- 
get qubit bias et , for a fixed control qubit bias Ec = and inter-qubit 
distance R — 300 nm. Locations of zero effective exchange are 
marked by Gray circles. Inter-dot bias is normalized to QD confine- 
ment and detuning is measured from the singlet anticrossing point. 



the /3 couplings is predominantly controlled by its respective 
qubit bias (as seen by the red dotted line in Fig.lTJb)). 

Complementing the target and control single-qubit Hamil- 
tonians with the interaction terms, Eq. ([T]i, we find that at cer- 
tain bias points exchange is turned off (gray circles in Fig.[T]5). 
We define a sweet spot in the 2D bias space {et, ec) as a point 
that satisfies two conditions: Jti^t, fc) — Jt ~ 2/3f — 0, and 
Jci^ti Ec) — Jc — 2/3c = 0. The first condition facilitates an 
idle position of the target qubit, or a pure X -rotation in the 
presence of a magnetic field gradient, 5ht. Similarly, the sec- 
ond condition is set to turn off the control qubit exchange, 
resulting in minimal interference with the target operation. 
For R < 1200 nm we find that for each qubit there are al- 
ways two such intersections, thus in the 2D bias space there 
are four possible working positions: two symmetric low- and 
high-bias points, and two asymmetric points. As explained 
below, the symmetric low- and high-bias points are optimal 
for single- and two-qubit operations, respectively. Fig. |2|a) 
depicts the dependence of the symmetric low-bias (solid red 
line) and high-bias (solid blue line) sweet-spot positions on 
the inter-qubit distance R. The corresponding entangling cou- 
pling, 7, at these bias positions is shown in Fig.|2b). For al- 
most the entire range, the low-bias 7 values exhibit a counter- 
intuitive increase with R due to the increased sweet-spot bias 
as R grows. The strong bias dependence of 7 overtakes its 
expected decrease with R for the considered range. 

Idle position and single-qubit gates. While the 7 coupling 
cannot induce relaxation of the two-electron qubit states, it 
generates dephasing of a superposition S-Tq state, and gate 
errors. Inspecting the two-qubit Hamiltonian, we find that 
without magnetic field gradients, 5ht = She = 0, when the 
two qubits are biased at a sweet-spot, preparing one of them 
in the state (|5) + \Tq))/\/2 results in a time evolution of 
its off-diagonal density matrix element: psTo (t) — 5 cos 27t, 
irrespective of the state of the other. Choosing the low-bias 
sweet-spot and R — 300 nm, we have 7 — 0.013 neV, corre- 
sponding to dephasing time of ~ 40/is. 

Since nuclear state preparation that generates a stable mag- 
netic field gradient requires pumping times between 60 ms - 
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FIG. 2: (color online) (a) Normalized target and control qubit bias 
at symmetric low-bias (solid red line) and high-bias (solid blue line) 
points vs. R. The dashed blue line depicts the high-bias position 
for the 7-corrected sweet-spot. The dotted black line marks the 
S(l, 1) — ^(O, 2) anticrossing bias point; (b) 7 coupling values at 
low-bias (red line) and high-bias (blue line) sweet spots. 



1 s 1 9], it may well be necessary to operate under a fixed Sh 
in both qubits. As long as She >> 7, the oscillation am- 
pUtude of the off diagonal element of the target qubit den- 
sity matrix reads 2{'y /6hc)^. Taking a conservative value of 
She — l/^eV, corresponding to a field gradient of ~ 40 mT, 
we find that at i? = 300 nm, target qubit errors are below 
10^^. Working with a finite Sht = 6hc = Sh, the sweet- 
spot is no longer an idle position for the target qubit, and 
the gate time for a 7r-flip X-rotation is 2 ns, with a lead- 
ing gate error of {n^/iShY ^ 4 x 10~^°, for Sh = l^eV 
1 1411 . If the Jt — point is used as an idle position, X- 
rotations can be erased either by waiting integer number of 
full periods, r — 2Trn/Sht, or by applying the pulse sequence: 
Uz{TT)Ux{oi)Uz{T!'), where a = ShtT is the acquired X- 
rotation during the wait time r. 

To generate Z-rotations with a finite Sh one can positively 
bias the target qubit so that J* >> Sh produces a nearly 
perfect Z-rotation. Note, however, that to maintain short 
gate times, Sh cannot be too small, thus a large excursion 
from the sweet-spot is necessary to obtain sufficiently large 
Jt- At this position, near the singlet anticrossing, the 7 cou- 
pling increases dramatically, generating substantially larger 
gate errors. In addition, a positively biased qubit may de- 
cohere faster due to environmental charge fluctuations. To 
stay near the Jt — point, one can perform a three pulse 
sequence: Uz{a) = Ue{x)Ux{(i>)Ue{x)r where ([/^(x) is a 
rotation about the 6'-tilted axis, = tan~^ ((5/t/ Jt), and the 
rotation angles x and are fixed by a and 9 1121 Il4ll . A Z- 
rotation by any angle a can be performed using this cycle, 
as long as \Jt\ > Sh i\9\ < 7r/4). Choosing 6 = 7r/4 to 
stay as close as possible to the Jt = point, we complete 
a three-pulse Z-rotation of a — tt in 5 ns, with a gate error 
smaller than 10~^. The dominant contribution to the Z-gate 
error reads: 2{'jg / Sh)"^ {3 — cos V2n), where jg is evaluated at 



the 0-tilted bias position. Although 70 increases with Sh due 
to the necessary larger bias sweep to the Jt — Sh position, 
the ratio 'jg/Sh is monotonically decreasing in the considered 
range, keeping the gate error very small for 6h < lOiieV il4ll . 

Two-qubit gates. We now consider the high bias symmetric 
point, where the qubit coupling 7 is large: 7 ^ Sh. Without 
magnetic field gradients, the two-qubit Hamiltonian is diag- 
onal in the computational basis state {SS, TqS, STq, TqTq}, 
and reads at the sweet-spot: H ~ diag(— 7,7,7, — 7). A 
CPHASE gate is conveniently obtained by tuning the bias in 
both double dots to satisfy Jq = 2(/3g + 7), q — t,c, at 
which H = diag(7,7,7, — 37). Letting the system evolve 
for time r^p — iih/A'^, the resulting gate, up to a common 
phase, is: diag(l, 1, 1,-1) = CPHASE. The position of the 
7-corrected high-bias symmetric point is shown by the dashed 
line in Fig.|2ta). 

Working with finite field gradients in both qubits, we find 
three different gate error scalings: (i) population errors with 
control in \Tq) ~ 5h/2^, (ii) phase errors ~ {5h/2^Y, and 
(iii) population errors with control in \S) ^ {5h/2^Y, the 
latter being negligible in the considered 5h range. We have 
performed numerical simulations of the time evolution of the 
two qubits under a controlled 7r-pulse, where the target qubit 
is initially prepared in the state: ( | S*) + 1 Tq) ) /\/2, and the con- 
trol qubit is in either \S) (target remains in its initial state) or 
|To) (target evolves to (IS*) - \To))/^). Fig.^a) shows the 
resulting CPHASE gate errors for R = 300 nm (7 = 20/xeF), 
which follow closely the above scalings up to 5h < IQ^eV. 
The small difference in the two phase errors with different 
control qubit states is due to the Hamiltonian structure at the 
7-corrected working point. 

To construct a CNOT gate, a Hadamard gate, Ur = 
{Ux{t^) + Uz{t^))/V2, should be applied to the target qubit 
before and after the controlled operation [15]. Uh can be 
implemented by biasing the target qubit to a position adja- 
cent to the low-bias sweet spot, where Jt = 5ht, for a dura- 



u 



tion TH — irh/ J Jt + Shf. We note that the control qubit 
should not rotate prior to the CPHASE gate. While it may 
be possible to employ a pulse sequence to correct such spuri- 
ous control rotation, here we take She = during the first 
Hadamard target gate. This may be implemented by, e.g., 
delaying the nuclear state preparation in the control qubit 
by th- This pulse sequence was simulated by discretizing 
bias and time steps, so that the actual switching times be- 
tween different bias points are taken into account. Fig. Ob) 
shows the resulting CNOT gate errors vs. Sht for R — 300 
nm. Population and phase errors interchange with respect 
to their role in the CPHASE gate, therefore the phase error 
with control in \To) is the largest here, scaling linearly with 
Sht/2'^. The results demonstrate increased errors for both 
small Sht (during the Hadamard gates), and large Sht (dur- 
ing the CPHASE gate). A minimum phase error of 4 x 10^^ 
is obtained for Sht = 0.017/iey, but with such a small gra- 
dient, single qubit gate times are on the order of 100 ns (see 
Fig-Elc)). Phase errors smaller than lO^'^ are obtained with 
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FIG. 3: (color online) (a) CPHASE gate errors vs. &h, for a target 
qubit initially at (jS) - |^o))/^/2, and i? = 300 nm. The solid blue 
(dashed green) line depicts phase error ~ (5/1/27) for control in 
singlet (triplet) state, and the dotted red line shows population errors 
^ 5/1/27, when control is in triplet. Population errors when control 
is in singlet are very small (not shown); (b) CNOT gate errors vs. 
hht , for an initial target singlet state (same results are obtained for an 
initial target triplet state). The solid blue (dashed green) line depicts 
population errors for control in singlet (triplet) state, and the dotted 
red line shows phase error when control is in triplet. Phase error 
when control is in singlet is very small (not shown); (c) CNOT gate 
time vs. Sht. Note that Tcp = 0.026 ns is fixed by 7, and most of the 
CNOT gate time is spent on the Hadamard target gates. 



0.004 < Sht < 0.07/iey. Other errors, and in particular sin- 
glet return errors, are considerably smaller. 

We stress that the opposing requirements for small single- 
and two-qubit gate errors, while operating under fixed Sh, are 
met by utilizing the high tunabiUty of the 7 coupling, whose 
values show more than six orders-of-magnitude difference be- 
tween the low- and high-bias sweet-spots for R — 300 nm. 
Taking R = 200 nm will increase the 7 ratio by another order- 
of-magnitude, resulting in a minimum phase error of 6 x 10~^ 
at Sht = Q.02fieV, and a wider range for which errors are be- 
low 10^'^ up to Sht ^ 0.65/iey. Reducing the inter-qubit dis- 
tance may, however, induce errors due to electron exchange 
coupling between the two double dots that are not included in 
the current analysis. Further optimization of gate errors and 
times is likely to be obtained by engineering the dot size and 
tunnel coupling within each double dot 11211 . 

Qubit dephasing due to fluctuating charge environment. 
Exchange-coupled spin qubits are vulnerable to dephasing in- 
duced by charge noise, due to the different c harg e distribu- 
tions of the singlet and triplet states 0131 lla llTll . The for- 
mulation we employed to find inter-qubit couplings can be 
used to quantify couplings of each of the qubits to nearby 
two-level charge fluctuators (TLFs) and evaluate the result- 
ing qubit dephasing. We identify two types of TLFs: (i) /3- 
coupled, where the charge jumps in and out of the trap, and 
(ii) 7-coupled, where the charge fluctuates between two sites 
in the trap. While /3-coupled TLFs have stronger coupling to 
the spin qubits, they are less abundant and locate farther from 



the QDs, near the 2DEG layer or quantum point contacts. 

Starting with a single TLF, we represent the qubit state with 
a Bloch vector, and solve for its precession under the pseudo- 
field: B(t) = {Sh,0,J -v^{t)) UM- ^(t) = ±1 represents 
a Poisson switching process, and v = (3,2j is the qubit-TLF 
coupling strength for /3-type and 7-type TLFs, respectively. 
Averaging over the stochastic process, the qubit signal de- 
phasing is given by: W{t) = (e"''*), with ip = v L (,{t')dt'. 
For free induction decay (FID), we find 1 19]: 



w^™w = ^E±f(i±A^)' + (f) 



,±rtit 



(3) 



where ^^ ^ 1 - (v/T)^ + 2iv/Ttanh{AE/2kBT), T = 
i(r+ + r_) is the average switching rate, r+/r_ = 
g-AB/feflT^ and AE is the TLF level splitting. To calculate 
signal decay during a general dynamical decoupling pulse se- 
quence we write ip — v J^ C(t,t')^(t')dt\ where C{t^t') is a 
filter function in the time domain flu]. For the spin echo (SB) 
protocol, C{t, t') = e{t')Q{t/2 - t') - e{t' - t/2)e{t - t'), 
where 0(i) is the step function, and the signal decay is found 
to be: 



-Tt 



WsB{t) = 
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(Ml. 



(a* 
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± 



(4) 



where fij^ (^j) is the real (imaginary) part of /i. 

In order to estimate two-spin qubit dephasing due to an en- 
semble of TLFs, we use W{t) = nr=fi wf{t) ]X^2^ W^{t), 
where n^ (n^) is the number of /3-type (7-type) traps. We 
estimate the trap density as ^ 120firn^^, and the two double- 
dot device area as ^ 0.75/im^, resulting in an active number 
of traps: n^ = 3, n^ = 9 [i21i] . In this estimate we assume that 
an area of radius ^ 150 nm around each qubit is depleted, so 
that traps in this area are not charge-active, setting up the max- 
imum qubit-TLF coupling strength. Trap parameters include 
level splitting of 10/ieF, center radius of 5 nm, and intersite 
distance of 20 nm, chosen to characterize (5-doped dopants in 
the insulator Finalh^ trap switching rates are taken between 

1ms < r-i < isjiai. 

We have used Eqs. ^ and (|4]i to calculate dephasing of a 
qubit at the low -bias sweet spot, by averaging over 10,000 sets 
of TLFs with random locations, orientations, and switching 
rates, within the parameter ranges specified above. In the FID 
case we obtain dephasing time of Tj* — 23 ns, comparable to 
measured coherence times (|5|]. For the SE protocol we find 
rf^ — 162/is, two-orders-of-magnitude longer than mea- 
sured SE coherence times |7], attributed to nuclear spin fluc- 
tuations, mediated by the hyperfine interaction ll20ll23ll . In the 
quasi-static regime considered here, FID dephasing times are 
dominated by the strongest fluctuator(s), and are thus mostly 
sensitive to qubit bias and qubit-TLS distance. SE dephasing 
times are dominated by the fastest fluctuator(s), and are sensi- 
tive to the maximal switching rate 1.24.] . We stress that for such 



a mesoscopic system, these results are not self averaging and 
may vary considerably between samples. For SE experiments 
with larger devices, sample-to-sample variability is reduced 
and dephasing times become comparable with those induced 
by the nuclei 1 14], thus both charge and nuclear environments 
should be considered in future device design and operation. 

In conclusion, we have quantified the couplings between 
two S — To qubits, identifying optimal working points. Our 
proposed scheme relies solely on interdot bias tuning at each 
qubit, and offers several advantages, including high fidelity 
single- and two-qubit gates under fixed field gradients, and 
idle positions at which the qubits are effectively decoupled. 
In addition, we have used our formulation to quantify qubit 
dephasing due to charge fluctuations. 

The author thanks Xuedong Hu and Hendrik Bluhm for 
helpful discussions and acknowledges funding from Research 
Corporation. 
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Supplemental material for: Electrically controlled quantum 
gates for two-spin qubits in two double quantum dots 



Coulomb couplings between two double-dots 

Here we outline the calculation of the Coulomb couplings 
between the two-electron orbital states associated with each 
qubit. Since all the gate operations considered in this work 
are performed by tuning the bias between the two QDs, the 
Hund-Mulliken orbital Hamiltonian must account for biased 
dot configuration. More details are given in [13]. 

We start by approximating the orbitals for the biased 
double-dot by those of two harmonic wells centered at (±a + 
eaB)x, where a is the interdot half separation, ub is the QD 
Bohr radius, and e = eEaB/fv^o is a dimensionless parameter 
corresponding to the interdot bias (normalized to the QD con- 
finement). We denote the orthonormalized (predominantly) 
left and right single particle orbitals by ^p±a, and use them 
to construct the two-electron spin basis states. These are 
the two doubly occupied singlet states, 5(2, 0) = ip-aip-a, 
5(0, 2) = ipaipa, and the separated singlet and unpolarized 

triplet states, 5(1, 1) = {ip-a^a + 'lpa'lp-a)/V2, To{l, 1) ^ 

{tp~ai^a — 4'aip-a)/\/^, where(i, j) indicates the number of 
electrons in each dot. Since the orbital Hamiltonian does not 
connect the triplet state with any of the singlet states, we only 
need to diagonalize the 3x3 singlet block, resulting in a hy- 
bridized singlet state. 

We note that both the polarized and doubly occupied triplet 
states are neglected in our analysis. Doubly occupied triplet 
states require one of the electrons to reside in an excited or- 
bital, thus their energy is far above the doubly occupied sin- 

n 
glets f5]. Utilizing Overhauser fields to generate a magnetic 

field gradient will likely create a small perpendicular compo- 
nent that may mix the Sz ~ computational subspace with 
the Sz — ±1 subspace , spanned by the polarized triplet states. 
With dh — IfieV, the transverse inhomogeneous field should 
not exceed Sh± = 50 neV for which leakage errors are found 
to be below 10~^ 11211 . thus polarized triplets can be safely 
neglected as well. 

For the Calculation of the inter-qubit Coulomb couplings, 
it is convenient to work in the orthonormal basis of symmet- 
ric and antisymmetric orbital combinations: i/ii — {tp-a + 
ipa)/V2, "02 = (V'-a — '0a)/\/2- In this basis, the Coulomb 
interaction operator between the two qubits is written as 



fijki = / drdr 



,p\A^)Pli{^') 



e r - 



(SI) 



V 2 (/nil + /2211 



in terms of the fijki operators: 

^TqTo = /nil + /1I22 + ./22II + /2222 
^ToS = "^^ToTo+A/'c (ai+a2 

— /1122 — /2222) + («! — 0-2) (/1II2 + /2212 
+ ./112I + /2221)] 

VsT, = Vt.To +^t {a\ + 4) [n/2 (/nil + /1122 

— /22II — /2222) + (fti — 0.2) (/12II + /1222 

+ ./2111 + /2122)] (S2) 

Vss = VtoT, + VtoS + VsT, + M^N^ {a\ + a* ) 
X {al + al) [2 (/nil + /2222 - /1122 - /2211) 
+ V 2 [a^ — flj) (/1211 + /2111 ~ /1222 — /2122) 
+ V 2 {al — flj) (/1112 + /1121 — /2212 — /2221) 
+ (oi - 4) (oi - 0-2) 
X (/1212 + /2121 + /1221 + /2121)] , 

where the superscript t (c) denotes target (control) qubit, 
fli, 02 are the 5(2, 0), 5(0, 2) components of the lowest lying 
hybridized singlet, respectively, and M is the normalization. 

The Coulomb interaction terms fijki are calculated by eval- 
uating the electrostatic energy associated with placing the 
control qubit charge distribution in the potential $* , that 
is due to the target qubit charge distribution. Denoting 
the charge, dipole, and quadrupole electric moments associ- 
ated with the qubit charge distribution by qij,Pij,Qij, re- 
spectively, we obtain fijki up to and including quadrupole- 
quadrupole order: 
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In deriving Eq. (IS3b we assumed both qubits are aligned with 
the a;-axis. The qubit dipole moments are thus in the x di- 
rection for both qubits, P^j = / drrpij{r)x. The quadrupole 
moments, Qij^^n) = J dr(3r™r„ - r'^5mm)pij (r), are found 
to have only diagonal elements {m — n). 

Explicit expressions for the various multipole moments are 
given in il3n . and are used to obtain the inter-qubit cou- 
pling terms, to quadrupole-quadrupole order with the follow- 
ing nonvanishing contributions: 



where Pij{r) = etp*{r)^j{r) is the electron charge density 
operator for either qubit, and i,j G {1,2} denote the qubit 
orbital state (symmetric or antisymmetric combination). The 
Coulomb matrix elements in Eq. dU in the main text are found 
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The first (second) superscript in each term denotes contribu- 
tion from the particular multipole moment: monopole (1), 
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FIG. SI: (color online) Multipole expansion contributions to inter- 
qubit Coulomb coupling vs. inter-qubit distance. Both double dots 
are biased at the 5(1, 1)-S'(0, 2) anticrossing point, (a) /? terms; (b) 
7 terms. 



dipole (2), and quadrupole (4) of the target (control) qubit 



charge distribution. The effective exchange contribution for 
the control qubit, /3c, is found from Eq. (IS4l i by interchang- 
ing the left and right superscripts. Symmetry considerations 
suggest that for identical qubits, equally biased, one has: 
/3j-' = /3;?', and 7'-' = 7-'*. Note that the dipole moment of 
a non-biased double dot vanishes, therefore there are no con- 
tributions to /3t that involve a non-biased target qubit dipole 
moment. Similarly (3c will have no control qubit dipole contri- 
butions for non-biased control qubit, and 7 will have no dipole 
contributions from either qubits if both are non-biased. 

Figure [ST] shows the dependence of the various multipole 
terms on the inter-qubit distance for the case of two double- 
dots equally biased at the singlet anticrossing point. The 
figure demonstrates the convergence of the multipole expan- 
sion in the considered distance range. We conclude that for 
R ^ 300 nm, used throughout most of this work, both j3 and 
7 couplings are well approximated by their respective leading 
contributions. 



Inter-qubit dynamics 

The two-qubit Hamiltonian, H — Hint+'Ht+'Hc, is^ written 
in the computation basis states: {SS, TqS, STq, TqTq}, where 
the left (right) script denotes the target (control) qubit state, as: 
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(S6) 



At the sweet- spot one has: 2/3^ = Jq, q = t,c. We consider 
dephasing of the target qubit, initially prepared in the state: 
{\S) + |To))/-\/2, by solving for the time evolution of the sys- 
tem of two qubits under the Hamiltonian, Eq. (IS6I 1, and tracing 
out the control qubit states. Notice that the term dephasing is 
used here and in the main text to quantify different time evo- 
lutions of S and Tq constituents in a superposition state due to 
inter-qubit coupling. This is not pure dephasing in the sense 
of loss of coherence due to entanglement with a bath. For a 
control field gradient satisfying 6hc ^ 7, which is easily met 
at the low-bias sweet-spot with a conservative magnetic field 
gradient of 40 mT, we find: 
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While oscillations are much more rapid (^ She) than in the 
She = case (~ 27), their amplitude is very small, thus the 
qubits are effectively decoupled. 



To gain better insight to the qubit decoupling bias points 
and their robustness against bias errors we plot in Figs.|S2ja,b) 
target qubit dephasing times vs. the control qubit bias, mea- 
sured from the Jc — 2(3c, low-bias position. We define de- 
phasing time as the time it takes p^^ (i) to drop to 50% of its 
original value of 1/2, due to coupling with the control qubit. 
In general, we note that the effects of the control qubit on 
target dephasing are minimal at a relatively wide control bias 
region around the sweet spot. This should be contrasted with a 
much higher sensitivity to bias errors at the high-bias position. 

When 6hc = 0, target dephasing at the sweet spot is iden- 
tical for both control states, but Fig. |S2f a) shows that slight 
detuning of the control bias results in complete qubit decou- 
pling when the control is in either singlet or triplet state. To 
explain this behavior we notice that at bias positions satisfy- 



ing Jq — 2(/3q ± 7) = in both qubits, Eq. (IS6b leads to: 
exp{~iHt) = diag(l, 1, e^^''''*, e^'''*) for the plus sign, and 
cxp{—iHt) = diag(e~^*^*, e^*^*, 1, 1) for the minus sign. As 
a result, the qubits are decoupled when the control is in a sin- 
glet (triplet) state for the plus (minus) sign, and the target qubit 
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FIG. S2: (color online) Dephasing times of target qubit vs. control 
qubit bias measured from the low-bias sweet spot, at i? = 300 nm. 
The control qubit is prepared in a singlet (blue squares) or triplet (red 
plus symbols) state, with magnetic field gradient of: (a,c) 5hc = 0, 
(b,d) She — IfieV. In figures (c) and (d) the sweet spot is shifted, 
satisfying Jq — 2{l3q + 7) = inboth qubits. 



Constructing .Z-rotations with finite 5h 

Here we use a variant of the bias sequence generating X- 
rotations that was proposed in 11211 . Our cycle consists of two 
working points: (i) the low -bias symmetric sweet-spot Jj = 0, 
where rotations are about the X axis, and (ii) a bias position 
slightly above the Jt = point, where Jt > 6ht and the 
rotation is about a 6'-tilted axis where 9 = ta,n^^{Sht/Jt). 
When \9\ < 7r/4, a Z-rotation by an arbitrary angle a can be 
generated by the cycle: 



U,{a)^Ue{x)UAme{x), 



(S9) 



where Ue{x) ^nd Ux{(j>) are rotation matrices about the 
tilted and X axes, respectively, given by: 
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The angles x ^nd are functions of 6 and a, and are found to 
be similar (but not identical) to those given in 1 12]: 



dephases as p^gj, [t) — i cos(47t) when the control is in the 
other state. Since /3c grows faster than Jc with €c, the inter- 
section Jc = 2 (/3c + 7) occurs at a negative bias, and simi- 
larly Jc — 2 (/3c — 7) occurs at a positive bias relative to the 
sweet spot position (see figure[TJb)). These considerations are 
numerically confirmed by correcting the sweet spot to satisfy 
Jq — 2(/3g + 7) = for (7 = t, c. Fig. [S2l c) shows qubit 
decoupling at this 7-corrected sweet spot position for singlet 
control state. When 5hc ^ 7, the qubits decouple at the origi- 
nal sweet spot (with errors of order ^ (j/Shc)'^) according to 
Eq. iSli . for both control states, thus shifting the sweet spot 
bias results in a similar shift in the decoupling bias position, 
as seen in Fig. lS2f d). 

Finally, eiTors in X-rotations performed at the Jt = bias 
position can be evaluated, using Eq. (IS6I ). Considering a tt- 
flip on a target qubit, initially prepared in a singlet state, we 
take T = w/Sht and find the leading contribution in a j/Sh 
expansion. At the end of the pulse, the gate error, given by the 
singlet probability is found to be: 
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2"f cos (^Shc/Sht) 
5ht{l-5hl/5h1) 



(S8) 



irrespective of the control qubit state. When Sht — She = 5h 
Eq. ( IS8b reduces to Ps{t) = {TT'j/2dh)'^, given in the main 
text. Notice, that when She = nSht, where n is an odd integer 
larger than 1, the leading contribution, Eq. ( IS8b . vanishes and 
we get X-gate error of the order {■y/Sht)'^. 
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where the plus (minus) sign in Eq. JS12I ) corresponds to pos- 
itive (negative) 9. In particular, a 7r-rotation about the Z axis 
is obtained with: 



X = arccos(— tan^ 0); (/> = ±2 arctan 



sin0 
Vcos 29 ' 



(S13) 



We consider a 7r-rotation about the Z axis by taking the sec- 
ond working point to satisfy 9 = n/A so that the bias sweep 
from the low-bias sweet-spot is minimized. The pulse dura- 
tions are set to: tq = TT/^/2Sht, and t^ = n/Sht, at the 
^-tilted point and low-bias sweet-spot, respectively. To sim- 
ulate this cycle we discretize bias and time steps, so that the 
actual switching time between the working points is taken into 
account. Bias switching should be nonadiabatic with respect 
to nuclear mixing, but slow as compared with the tunnel split- 
ting of the hybridized singlet states. Bias sweep rates < 2 
mV/ns, comparable to experimental values ||5|], do not con- 
tribute appreciably to gate errors presented in this work, and 
can be further corrected by pulse design. 

Fig. ISST a) depicts the target off-diagonal density matrix 
element during a vr-rotation about the Z axis from (jS*) + 
\To))/V2 -^ {\S)-\Ta))/V2, using the sequence in Eq. ^. 
Taking 6ht = She = Sh we use Eq. ( IS6b to find the devia- 
tion of the the target off diagonal density matrix element from 
— 1/2 at the end of the cycle. The leading error terms in a 
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FIG. S3: (color online) (a) tt -rotation about the Z axis using a three- 
pulse sequence. 5ht — l^eV and R — 300 nm. The first and third 
rotations are performed at a bias position satisfying Jt — She (9 — 
7r/4); (b) Gate error as a function of Sh, evaluated from numerical 
simulation (solid blue line) and from Eq. ( IS14b (dashed red line). 



"fx/5h, ^g/5h expansion, where ^^ ije) is evaluated at the 
sweet-spot (6'-tilted) position, are found to be: 
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and are compared with the numerical simulations in 
Fig. lS3r b). The proposed sequence, operating around the low- 
bias sweet-spot, is shown to generate high fidelity ^-rotations 
at a reasonably short gate time. 



Alternative CNOT gate 
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FIG. S4: (color online) Two-spin qubit dephasing due an ensemble 
of TLFs, for two device sizes. Figs, (a) and (b) show results for 12 
TLFs, and Figs, (c) and (d) show results for 1 12 TLFs. Two qubit ma- 
nipulation protocols are considered: free induction decay in Figs, (a) 
and (c), and spin echo in Figs, (b) and (d). In all figures the solid 
blue lines correspond to an average over 10,000 sets of TLFs with 
random locations, orientations, and switching rates, within the pa- 
rameter ranges specified in the main text. The dashed green and red 
lines mark the distribution widths of the signal decays. 



double dots in opposite directions, so that the configuration: 
S{0, 2) — 5(2, 0) is obtained. This bias position is very far 
from the low-bias symmetric point that is proposed in the cur- 
rent work for single qubit operations, and could be experimen- 
tally challenging. 



For demonstration purpose an effective CNOT gate for an 
initial target qubit in either \S) or |To) can be implemented by 
biasing both qubits to the point Jg — 2{/3q — 7) = 0, g = i, c, 
close to the high-bias symmetric sweet-spot discussed in the 
main text. At this bias, the diagonal part of the two-qubit 
hamiltonian is (—87,7,7,7), ^"d "f ^ Sh < IjieV. Thus, 
the target state will effectively rotate about the Z-axis when 
the control is in \S) and about the X-axis when the con- 
trol is in |To). Applying this pulse for t = n/Sh to a two- 
qubit system, with the target initially in either of the compu- 
tational states (but not in a superposition of the two), results 
in a CNOT gate, with gate time inversely proportional to 5h. 
The degeneracy in the Hamiltonian at the above bias point in- 
duces gate errors that scale linearly with 5h/^, and therefore 
are rather large for R = 300 nm (7 = 20iJ,eV). A better bias 
position is the slightly asymmetric point: Jt — 2(/3t — 7) — 0, 
Jc — 2/3c — 0, at which the diagonal part of the two-qubit 
hamiltonian is (—27, 27, 0, 0). At this position gate errors are 
found to scale as (Sh/j)'^, and are small over a wide range of 
dh values. 

CNOT gate can also be implemented using the scheme 
proposed in 11211 . but this scheme requires biasing the two 



Effect of TLF number on qubit dephasing 

Here we demonstrate the effect of the device size (and 
hence the number of charge traps) on qubit dephasing and 
sample-to-sample variability. Figs. |S4] (a) and (b) show cal- 
culated FID and SE dephasing, respectively, for a device of 
area O.Tbiim? (np — 3, n^ ~ 9), considered in the main text. 
The solid blue lines depict average decay over 10,000 random 
TLF sets, and the dashed green (red) lines correspond to max- 
imally (minimally) affecting sets. The results show a large 
sample-to-sample variability, in particular for the SE decay, 
which is governed by the wide switching rate distribution, 1 
ms< T-^ <1 s. Figs. [Si (c) and (d) show FID and SE de- 
phasing, respectively, for a device of area 6iim^ (np = 28, 
Uy — 84). The FID dephasing time, Tj* = 17 ns, and variabil- 
ity are hardly changed with respect to the smaller device. In 
contrast, the SE decay shows an order-of-magnitude reduction 
in dephasing time, r|^ = 2G^s, and much more reproducible 
decay, suggesting that for ntrap ^ 100 results become self 
averaging. 

We note that in both qubit manipulation protocols consid- 
ered here, the noise generated by the charge environment is 



9 

non-Gaussian, and therefore cannot be described in terms of The dependence of qubit dephasing on TLF parameter distri- 
its spectral function alone. For larger devices, where the num- butions, qubit characteristics, and magnetic field, and its scal- 
ber of fluctuators is of the order of 100, the distribution of ing with the number of qubits, as well as qubit dephasing un- 
coupling strengths and switching rates can be averaged, and der multiple-pulse dynamical decoupling sequences, will be 



the resulting noise is dominated by a large number of weakly presented elsewhere 124 1 
coupled TLFs, demonstrating characteristics of 1// spectrum. 



